Scattering in plane-wave backgrounds: infra-red effects and pole structure 
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We consider two aspects of scattering in strong plane wave backgrounds. First, we show that the 
infra-red divergences in elastic scattering depend on the structure of the background, but can be 
removed using the usual Bloch-Nordsieck approach. Second, we analyse the infinite series of shifted- 
mass-shell poles in the particle (Volkov) propagator using lightfront quantisation. The complete 
series of poles is shown to describe a single, on-shell, propagating particle. 



I. INTRODUCTION 

Strong external fields affect many aspects of gauge 
theories. In QCD, magnetic fields affect the vacuum 
PQ, phase diagram [2 j, electric dipole moments [3 and 
the quark-gluon plasma [4 j. In QED, scattering pro- 
cesses in the fields of intense lasers currently attract quite 
some interest [5j [6] , with the aim of investigating both 
nonperturbative effects [7] and beyond-standard-model 
physics [8[ [9]. Modelling the laser as a plane wave al- 
lows scattering amplitudes to be calculated for arbitrarily 
strong fields because the fermion propagator in a plane 
wave is known exactly; this is the Volkov propagator [TO] . 

In this paper we will consider the propagation of quan- 
tum particles in strong plane waves. While the basic re- 
sults for scattering in plane waves were given in the 1960s 
[TTHT3] , those calculations assumed either monochro- 
matic waves or crossed fields (constant plane waves). 
Both of these fields are of infinite extent and are there- 
fore rather special cases; two statements related to them 
will be examined below. 

We begin by considering the infra-red ( "IR" ) structure 
of processes in strong plane waves, focussing on soft cor- 
rections to elastic scattering. In a crossed field, the dif- 
ferential probability of photon emission scales like 1/uj 2 ^ 3 
rather than as 1/uj as in bremsstrahlung. Thus the log- 
arithmic divergence of the IR becomes an integrable sin- 
gularity. That this happens in a field which never van- 
ishes is potentially interesting because IR problems origi- 
nate in the (incorrect) assumption that the QED coupling 
switches off at large distances [14] [15] . This weakening of 
the divergence actually comes at the expense of admit- 
ting unphysical large distance behaviour, see [TO] , but it 
has raised the question of whether a partially nonper- 
turbative treatment of plane wave backgrounds can offer 
insight into the IR problem [T7HT9] , which is is still un- 
der active investigation [20-25 . It is also possible that 
background fields lead to new problems; it has been sug- 
gested, for example, that IR divergences do not factorise 
in pair-creating backgrounds [26] . 

In the second part of the paper, we turn to the basic 
building block of elastic scattering, the electron propa- 



gator. Our focus is on the poles of the propagator, par- 
ticularly in monochromatic fields. These are still used as 
an intuitive basis for more general calculations in finite 
pulses, and this has led to some debate concerning the 
'intensity dependent mass shift' [27l [28] and [29] [30]. The 
mass shift (or rather, its effect) can be seen in the spec- 
trum of undulator radiation [3T] , but its theoretical def- 
inition is more elusive. When it appears, the mass shift 
leads to poles in the propagator away from p 2 = m 2 , sug- 
gesting the presence of heavy states. We investigate this 
by directly constructing the quantum states of a particle 
in a plane wave and also by resumming the pole contri- 
butions to scattering amplitudes. 

The paper is organised as follows. After a brief review 
of previous results, we discuss in Sect, [il] loop and soft 
emission corrections to elastic scattering in plane wave 
backgrounds, and the canc ellation of IR divergences fol- 
lowing [32, 33 . In Sect. Ill we discuss the propagator and 
the quantum states using lightfront quantisation, and ex- 
plicitly relate the poles of the propagator to the ordinary 
mass shell condition. We conclude in Sect. |IV| Appendix 
[A] collects useful results on the propagator and the nor- 
malisation of ^-matrix elements in external plane waves. 
Appendix [B] contains the details of our IR calculations. 



A. Conventions and review 

Consider a classical particle in a plane wave F^ v (fa) 
depending on <j) = k.x with k 2 = 0, lightlike. We take 
k.x = ux + , lightfront time. (Recall that x ± = x° ± 
x 3 , x ± = {x 1 ,^ 2 }, and x ± = 2x T .) We consider finite 
duration fields, for which F^ v vanishes before some fa and 
after some <j)f. The lightfront structure of the plane wave 
means that all initially present particles enter (leave) the 
wave at the same lightfront time fa ((f) f). A particle 
entering the wave with momentum p^ has a subsequent 
kinematic momentum 7r M given by 



7i>0; fa)=P^~ eC^ifa + k f 



2ep.C(fa-e 2 C 2 (fa 
2k.p 



, (1) 



in which C^(fa) is the integral of the electric field strength 
from the initial to the elapsed lightfront time: 
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and C± = 0. No gauge potential is employed in this 
derivation, see [34]. We now set fa = 0. We are in- 
terested in 'unipolar' pulses for which the integral over 
the entire electric field is nonzero [35], and which can 
be taken as crude models for fields which provide vac- 
uum acceleration [36-38 . This means that after F^ v has 
switched off, C M becomes constant, and we can identify 
it with its value at <\> = oo, so we write C^{fa) = for 
4> > (j)f. Using (fil), we then see that the difference in mo- 
mentum for an electron passing through the wave obeys 
(tt(p; oo) — p) 2 — e 2 C^ < , which has the right sign for 
scattering. When Coo = 0, there is no net acceleration 
and ir(p; oo) = p. Note that C^ is equal to the Fourier 
zero mode of the field strength [T6] . 

Regarding gauge invariance in theories with plane wave 
backgrounds, there is nothing to discuss classically: par- 
ticle motion is described by the Lorentz, LAD or LL equa- 
tions which depend only on F^ v . Quantum mechanically, 
one must choose a potential. For plane waves, the almost 
universal choice is to take (up to a constant) 

^background = e(J ^ > (3) 

with C M from This choice, often made implicitly, 

makes the physics manifest. Now that the gauge is fixed, 
final probabilities will depend on C M . We will see this 
in the appearance of 71-00 := 7r(p; 00) in the probabili- 
ties below. Using ([2|, though, we can always rewrite 
such results in terms of F^ V1 which secures gauge invari- 
ance. (This explains the apparent result that quantum 
processes in plane waves "depend on the vector poten- 
tial characterising the laser field, not on the electric field 
component" [39].) 

II. THE INFRA-RED SECTOR 

We now consider the infra-red behaviour of scattering 
processes in plane- wave backgrounds, through the use of 
several examples. The details and derivation of all our 
IR results may be found in Appendix B. We work in the 
Furry picture, in which the background field is treated ex- 
actly (i.e. without recourse to perturbation theory). We 
begin by recalling that the soft contribution to the prob- 
ability of single photon emission, or 'nonlinear Compton 
scattering' see Fig. [I] and [40J, is [16] 

F-=-e 2 [^L±(l^_P_Y (4 ) 
• J (2ny2l U^oc I.V) ' K ) 

using dim reg in d > 3 dimensions, and Y is understood 
to carry an upper cutoff corresponding to, say, detector 
resolution. This is logarithmically divergent in d = 3 
when Coo 7^ 0, but vanishes when Coo = 0- Hence, the 
presence of this IR divergence depends on the properties 
of the background. 

IR divergences typically arise when virtual particles 
come close to the mass shell. The field-dependent IR di- 
vergences in plane- wave backgrounds arise as follows. At 
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FIG. 1. Nonlinear Compton scattering of a soft photon, 
at tree level. Double lines indicate the background-dressed 
Volkov propagator. 

each vertex, the structure of the background allows the 
x~ and x 1 - integrals to be performed immediately. The 
p + integral can be performed using the residue theorem, 
which restricts the remaining x + integral by introducing a 
lightfront time-ordering. One then sees that it is the large 
lightfront time parts of these integrals which yield singu- 
larities. In other words, our IR divergences essentially 
arise from the background- free regions of spacetime, be- 
fore and after the pulse. See Appendix B. 

A. Double Compton scattering: hard-soft 
factorisation 

Consider now the emission of two photons from an elec- 
tron in a plane wave, i.e. the process 

e-(p) in ^ eT e-(p')+ 7 (*')+7(0, (5) 

as recently investigated in [4TH43] . Assume that the pho- 
ton with momentum / M is soft. It can be emitted from 
either the incoming or outgoing leg, and the ^-matrix 
element takes the form 

I p p f \ 

Sfi = ee so ft. -j 7—7 )S NLC (p p',k') , (6) 

\l.p l.p'J 

in which Snlc is the ^-matrix element for nonlinear 
Compton. This is the expected form of a soft correc- 
tion to a hard scattering process; two-photon emission 
becomes degenerate with nonlinear Compton when one 
of the emitted photons is soft. The soft divergence im- 
plied by (|6| is independent of the structure of the back- 
ground. This is an example of a general result (see the 
appendix) : the structure of the plane wave has no impact 
on the hard-soft factorisation of IR divergences, or the 
severity of those divergences, which give the usual l/s 1R 
poles in 4 + 2e IR dimensions. (See [44, 45 for examples in 
crossed fields.) Essentially, taking the soft limit removes 
the background-field dressing from emission vertices, and 
factorisation proceeds as in QED without background. 

However, the soft divergence implied by |6| is not the 
highest order divergence in two-photon emission. Rather, 
this comes from the case in which both photons are soft; 
two-photon emission then becomes becomes degenerate 
with elastic scattering, see Fig. [2] The IR divergent part 
of the emission probability in this case is 

pglbf ddl fjEgL_P.yi 2 = ly2 (?) 
2 [ J (2tt)«2Z Woo l.p) \ 2 • y > 
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FIG. 2. The probability for elastic scattering. The grey dot 
denotes all loop corrections. 



Unlike in (J6|, this divergence does depend on the struc- 
ture of the background. For C^ ^ each integral con- 
tributes the same divergent term; the leading singularity 
is therefore l/e? R . 



B. Elastic scattering 



As shown in Appendix B (leading to (B21)), the elas- 



tic scattering probability including the soft contribution 
from all loop orders, see Fig. [2j is 



exp 
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d d l 1 



{2ir) d 2l V-7T, 



P_ 
l.p 



(8) 



which defines X. This expression has two parts. The 
exponential is the all-orders soft loop contribution which 
is log divergent in d = 3 when C^ ^ 0- This multiples a 
'1' which is, see the appendix, the exact tree- level prob- 
ability of elastic scattering. (That P = 1 at tree level is 
already a sign that something is wrong.) 

These problems are related and their resolution is clear: 
it is not possible to observe 'elastic scattering' alone, due 
to the potential emission of soft, unobservable, photons. 
When we calculate the inclusive probability, IR diver- 
gences coming from soft emission should cancel those 
coming from the loops. We turn to this now. 



C. IR cancellation 

Generalising (|7|), the tree level probability of emitting 
n soft photons is Y n /n\. The all orders soft-loop contri- 
bution to each of these processes is e x , as in ml What 
can be observed in an experiment is the sum of 1) the 
probability for elastic scattering and 2) the probabilities 
for emission of any number of undetected soft photons. 
This sum, see Fig. [3j is the measurable probability of ob- 
serving scattering of the electron without photon emis- 
sion, and its IR part is 



P(e- ^e~) = e x -1 



(r + |r 2 + ...) 



^ e X+Y 



(9) 



We see that the IR contributions factorise. Comparing 
Q and ([8]), we see also that X = — Y, so that 



FIG. 3. The measurable probability of 'scattering without 
emission' is the sum of the probabilities for elastic scattering, 
and the probabilities for the emission of arbitrary numbers of 
soft (unobserved) photons. 



and the leading field-dependent soft divergences cancel 
to all orders. Two remarks are now in order. 

First, we have followed [33] and considered only the 
divergent IR contributions to amplitudes, showing that 
these cancel. In order to extend our results to the com- 
plete amplitudes, i.e. in order to include the IR finite 
parts, one can instead follow the method of [32]. (See 
[44U47] for various loop calculations.) 

Second, there are no purely soft divergences in a 
crossed field [12]. In such a field, all particles are acceler- 
ated to the speed of light. Further, the literature results 
assume that the particle is also initially moving at the 
speed of light. This leads to the replacements {7r,p} — ^ k 
in ([7]) and ([8]), so that X and Y vanish individually. If 
a field which is nonzero and constant for a long but fi- 
nite time is used instead, the Coo dependent logarithmic 
divergence reemerges [16] . Thus, the 'improved' IR be- 
haviour in crossed fields comes at the expense of intro- 
ducing unphysical large distance behaviour. It does not 
seem to say anything about the large distance structure 
ofQED. 



D. Indistinguishable processes 

Suppose now that C^ = 0, i.e. there is no vacuum ac- 
celeration, and consider again nonlinear Compton scat- 
tering. The final state integrals are now IR finite, and 
one can integrate over all photon momenta to calculate 
the 'full probability'. This number, though, is not mea- 
surable. Even when processes are IR finite, one must still 
account for indistinguishable processes. 

Physically, this issue originates in the nonzero fre- 
quency range of a pulsed field. If we parameterise this 
frequency range as sou where s is real, then each s can 
produce photons with a frequency uj' s1 where (taking the 
incoming electron to be at rest for simplicity) [16] 



1 



1 



-(1 - cos (9) . 



(11) 



-) = 1, 



(10) 



Each uj' s is bounded below, u' s > ujs, but this lower bound 
can extend down to zero frequency in a pulse, even if 
there is no support at zero frequency itself, i.e. if C^ = 
[16] . Suppose, then, that an experiment can detect only 
photons of frequency higher than ujq. The measurable 
probability for one-photon emission is then the sum of 
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that for nonlinear Compton scattering with uj' > c^o, to- 
gether with the sum of all n-photon emission probabilities 
in which n — 1 photons are soft, u/ < ujq. 

For a discussion of indistinguishable processes and the 
structure of the Volkov propagator, see [48]. The essen- 
tial part of the amplitudes calculated above is in fact just 
this propagator: elastic elastic at tree level is given by the 
double amputation of the propagator, see Appendix [Aj 
and this is multiplied at higher orders by soft corrections. 
We therefore turn now to the propagator itself. 



III. THE POLES OF THE PROPAGATOR 

The Kallen-Lehmann representation of the two-point 
function provides a natural way to identify the mass of 
single-particle states, as the location of the poles in mo- 
mentum space [49] . One of the most discussed properties 
of the Volkov propagator is its infinite series of poles at 



(p-lk) 2 =mt 



(12) 



where m 2 = m 2 (l + Oq) is the shifted mass of a particle 
in a plane wave of intensity ao [34, 50 , 51 . There is no 
pole at the ordinary mass shell. The I = pole, p 2 = 
might suggest a change in the rest mass, while the other 
poles describe such a heavy particle absorbing 'multiple 
photons' from the background. The literature contains 
many interpretations of the poles. They were originally 
considered to be problematic [52 , but that higher loop 
corrections would regulate the poles by giving them a 
finite height and width [53]. The fact that the poles are 
discrete has lead to the claim that the spectrum of a 
particle in a plane wave is discrete [52j [54] . See [55] for 
an interpretation of the mass-shift m* as a finite mass 
renormalisation, and [56 for an interpretation in terms 
of effective potentials. 

Our focus is on this pole structure. Spin effects do not 
impact our discussion so we restrict to a scalar particle. 
Propagator poles are usually tied to the optical theorem 
and intermediate states, so we will begin by constructing 
the quantum states of a particle in a plane wave explicitly. 
We will see that loop corrections are not, in fact, neces- 
sary, hence we turn off the QED interaction. We there- 
fore consider only a scalar particle in a plane wave; since 
this background cannot spontaneously produce pairs [57] , 
particle number is conserved and the theory is simple. 
However, the equal-time quantisation of this theory does 
not seem to appear in the literature. The principle diffi- 
culty lies in proving certain orthogonality relations of the 
Volkov solutions [58 , 59 . The reason for this difficulty is 
that the background singles out preferred lightlike direc- 
tions. It is therefore natural to tackle problems in plane 
waves using lightfront quantisation, as suggested in [60] . 
We will show here that the theory is trivially quantised on 
the lightfront. We will then recover the Volkov propaga- 
tor as a lightfront time-ordered product and investigate 
the poles. The explicit equal- time quantisation of the 
theory will be performed in the forthcoming article [48 . 



A. Particle states in plane waves. 

A complex scalar in a plane wave obeys the equations 
of motion (D 2 + m 2 )(p = 0. The general solution is 

, x f d 3 p 0(p_) . x 7 + . . . oX 

<PW = J ( 27r ) 3 2p a?(Pp ^ + P^ p W ' ( ' 

where p := {p±,P-}. The functions (p p are scalar Volkov 
solutions obeying the initial condition (p p (4> = 0) = 
e -zp.x ge ^. i s therefore the initial data at 

= 0, and the Volkov solutions recover the classical kine- 
matic momenta of a particle in a plane wave via 



iD^{(j))ip p {x) =7T At (p;0)^ p (x) . 



(14) 



To quantise, we impose the lightfront commutation rela- 
tion (LCR) EI], 



[tp(x),2d-tp\yj\\ x + 



i5^-{x-y). (15) 



Using (113]) to calculate the left hand side of (15), the 



complicated exponentials in the Volkov solutions cancel 
immediately, and the commutator reduces to that of the 
free lightfront theory. The LCR is obeyed if 



[op.oj] = [&P.&J] = (27r) 3 2p_^ 3 (p - q) 



(16) 



which are the usual free-field commutators on the light- 
front. The particle interpretation of our theory is as fol- 
lows. We define the lightfront vacuum | ) to be annihi- 
lated by the a p and b p as usual, and then the first excited 
states are 



|p):=aj|0), and |p):=6+|0), 



(17) 



which we will now see are one-particle/antiparticle states 
respectively, just as in the free theory. Using the en- 
ergy momentum tensor, we combine the normal ordered 
Hamiltonian and kinematic momenta into 



n M (0) = T_ M = J d 3 x d_^D^(x) 
'd 3 p0(p_ 



(18) 



(2yr) 3 2p_ 



where ir u 



The states 







then have time- 



dependent energies and momenta given by 

n M (0)b) = 7i>(p;0) \p) , 

n M (0)|p)=7f M (p;0) \p). 



(19) 



The excited states therefore carry the time-dependent 
momenta of the classical theory. So, and create 
particles which, on the initial surface = 0, have on- 
shell kinetic momenta p^. This labels the continuous 
spectrum of states. At all subsequent times, the states 
carry momentum 7r M (p; 0) with tt 2 = m 2 . 
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B. The propagator 



The well-known quasimomentum obeys q = m (1 



2^ 



yielding the shifted mass. The Volkov solu- 



We have seen that the one-particle states have mass tions in this field are 



m , and that this spectrum of states is continuous. The 
background-dependent structure in the field operators 
(essentially the Volkov solutions) simply describes the ac- 
tion of the Lorentz force on a particle in a background 
field. How should this be reconciled with (12) which 
might suggest that additional "heavy" states appear? To 
answer this we turn to the propagator, which is the light- 
front time-ordered product of the fields (13): 



G(x,y) = (0\TMx)<pHv)\0 
M 3 p%-) 



+ 0(y + - x + )(p_ p (x)(pi p (y) . 



(20) 



The lightfront time ordering may be made covariant by 
using an ie prescription to yield 



G(x,y) = i 



-ip.(x-y)-^ f 2eA.p-e 2 A 2 

d 4 p e k - y 



. f d 4 ; 



(2tt) 4 p 2 
d 4 p <p P (x)<pUy) 



(21) 



)4p 2 



This is the Volkov propagator [62j [63] . In the second 
line, we have written the Volkov solutions ip p (x) with a 
script label to indicate that p^ is arbitrary, i.e. p 2 ^ m 2 
in general. Performing the integral over p + (on which 
cp p depends trivially), the simple pole returns us to (20), 
putting the initial momentum and therefore the ki- 



netic momentum 71^, onto the mass-shell. 



Since we can write p 



m 



- m 2 



we see 



that ( [21] ) is a spectral representation of the explicitly 



time-dependent operator D 



rn 



It is not, in contrast 



to the free propagator, a Fourier representation, which 
is available only for simple fields such as monochromatic 
waves. Fortunately, this is just the case of interest. So, 
consider a a circularly polarised, monochromatic field 



CM 



I 1 sin ( 



1 + I COS ( 



(22) 



with e 2 l\P 



m l affi 3 , l J .k 



appears here and is equal to eE r 



0. The intensity ao 
s /muj. We define 



re — i 



: l\p 
k.p 



l\p 
k.p 



and 



Ik.p 



ip p (x) 



-iq.x — ir sin(</>— 0) — ir sin(^) 



(24) 



(The final term in the exponent follows from initial condi- 
tions, but is usually dropped in the literature. Our results 
hold in either case.) The Fourier transformed propagator 
can be constructed directly from (20) or (21) and coin- 
cides with the known result [55, 62[ I63j 



GmonQAp) = ^2 J n+i{r)Ji(r)e v 

n,iez 

(2tt) 4 £V -p-nk) 



(25) 



(p — Ik) 2 — m 2 + ie 



The delta-comb structure is due to the periodicity of the 
background [55] |64] . Since the Bessel functions are ev- 
erywhere regular, we see that when p^ and p'^ are such 
that the n th delta-function has support, we recover the 
infinite series of poles (12). 



C. Pole resummation 

Writing a Feynman amplitude in momentum space, in- 
ternal lines become the Fourier transform of the propaga- 
tor. From (25), we see that one then hits a "resonance" 



for particular values of I and n, each of which corresponds 
to taking particular values of energy-momentum from the 
backgrouncQ In general, all the poles terms contribute 
to a given amplitude, and there is no obvious way to 
single out a particular pole over any other. We therefore 
consider the sum of contributions from all the poles. This 
can be extracted by taking real and imaginary parts, i.e. 
applying the standard result 



1 V 

= — it:5(x) H 

x + ie x 



(26) 



to (25) and retaining only the delta function terms; call 



(23) this part of the propagator R. In position space, one has 



R(x,y) 



d A (r> n') 

\Jr'>y J ip.y-\-ip .x 

(2tt) 8 



E 



Jn+^J^ry^ {2n)H{p' -p-nk) nS((p - Ik) 2 - ml) . 



(27) 



Related structures are seen in the lower order three-point pro- 
cesses of strong-field QED; due to the periodicity of the back- 
ground, the emission rates take the form of a diffraction pattern. 



When the momentum transfer over a cycle is a multiple of the 
laser frequency, there is a peak in the emission rate which is 
analogous to a patch of constructive interference |60l I65| . 
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Performing the p' integrals, changing variables n — »• s = 

2 

n + Z and p — » p — lk + fc, and summing over the Bessel 
functions using 



J s (r)e- is ^ k - x -^ = exp ( - irsin(/cx - tp)) . (28) 



gives the final result 
d 4 p 



,,) = / 



n5(p 2 - m 2 )<p p (x)<p;(y) , (29) 



which is a sum over on-shell Volkov wavefunctions (24). 



The same result is obtained directly from (21 ) by sending 



p A 



16 



7r5(p 2 — rn 2 ) . 



(30) 



Hence, the total contribution of the infinite series of poles 
is to replace the propagator by an integral over all real, 
on-shell intermediate states. In other words, the poles 
contribute the 'imaginary parts' one obtains from cutting 
the propagator, in the sense of the optical theorem. This 
confirms that the off-shell poles do not describe heavy 
states. The reason that the usual Kallen-Lehmann inter- 
pretation does not go through directly is that its deriva- 
tion assumes Poincare invariance of the theory. This is 
explicitly broken by the presence of background fields 
(though covariance is not). Related to this, the Fourier 
transform cannot be interpreted in the same way as in 



a free theory since canonical momentum (the Fourier 
variable) and kinematic momentum are not equal. See 
[66] for analogous statements and an investigation of the 
Kallen-Lehmann representation in AdS space. 

Our treatment of the poles has been formal: the poles 
do play a role in the detailed structure of emission rates 
[42] , and they collectively describe the regime of momen- 
tum exchange in which sufficient energy is taken from the 
background to put normally virtual (intermediate) par- 
ticles onto the mass-shell. Associating individual poles 
to physical states is misleading, though. The confusion 
arises because in monochromatic (periodic) waves, the 
one-particle states are also eigenstates of cycle-averaged 
momentum operators, with constant eigenvalues equal 
to the quasi- moment a, which square to the shifted mass. 
One then speaks of the quasi-momenta as the 'good quan- 
tum numbers' of the system [67]. However, such nonlocal 
operators do not tell us much about the states. A particle 
in a background field represents a time-dependent prob- 
lem and so 'eigenvalues' are in general time-dependent, 
as in (19). 



The results above extend to general plane waves as 
follows. Given Volkov solutions ip p (x) in a particular 
plane wave, we define R(x,y) as in (29) and introduce 
the Fourier transform T via 



/ 



ds 



<p p (x) = e~^ x / ^ e""*-* r.(p) , (31) 
from which we obtain the implicit Fourier transform 



G(p',p) = 

It is then trivial to check that 
R(p',p) -- 



dl dn 



r„ +; (p)rr(p) (27r) 4 <5V -p-nk) 



(p — Ik) 2 — m 2 + ie 



dl dn 



T n+l (p) r,*(p) (2tt) 4 «5V -p-nk) Re 



(p — Ik) 2 — m 2 + ie 



(32) 



(33) 



The sum over on-shell intermediate states, R, is therefore given by taking the real part of the free propagator buried 
inside G. In this context the contributions from the background act as a 'dressing' which encodes the time dependence 
of the system, but the plane wave does not change the fundamental particle content of the theory. 

I 



IV. CONCLUSIONS 



Hard-soft factorisation of scattering processes in QED 
goes ahead in the presence of a plane wave background 
field of arbitrary strength and shape. The factorisation is 
not sensitive to the structure of the background. The im- 
plication is that, just as in ordinary QED, IR divergences 
in plane wave backgrounds exponentiate and cancel from 
measurable processes. A rough explanation of why is 
as follows. A scattering process in a fixed background 
obeying Maxwell's equations can always be rewritten as 
scattering between asymptotic, coherent, photon states. 
As these are free-theory states, the scattering process is 



equivalent to a sum over ordinary QED process with all 
numbers of photons. Hence, if the IR divergences cancel 
in QED, they should also cancel here. (Compare [26J, 
which suggests that pair-creating backgrounds not obey- 
ing Maxwell may lead to non-fact or isable divergences.) 

However, in those processes which are entirely soft, 
e.g. elastic scattering, the IR divergences depend on the 
Fourier zero mode of the field strength (if this is nonzero, 
the pulse can transfer net energy to a classical particle 
passing through it). Nevertheless, we have shown that 
the soft IR divergences in loop corrections to elastic scat- 
tering are cancelled by divergences coming from multiple 
soft photon emissions, as normal. 

The structure of the plane wave leads to the appear- 
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ance of light front time-ordering in scattering amplitudes, 
and IR divergences then arise from the large lightfront 
time regions before and after the pulse. This suggests 
that the natural setting for strong field QED is light- 
front quantisation, as employed in [60]. We used light- 
front quantisation to explicitly constructed the quantum 
states of a particle in a plane wave (which are contin- 
uous, not discrete as previously claimed), and recover 
the Volkov propagator as a lightfront-time-ordered prod- 
uct. We showed that the shifted mass-shell poles in the 
Fourier representation of the propagator actually corre- 
spond to going onto the ordinary mass-shell. The reason 
that the poles do not correspond to particle masses is due 
to the explicit breaking of Lorentz invariance induced by 
the background; the presence of the laser means that the 
Fourier variable does not coincide with physical momen- 
tum, and hence the Kallen-Lehmann interpretation of the 
poles does not apply to the Fourier representation of the 
propagator. It seems more natural, therefore, to talk of 
the mass-shift in terms of its observable effects, namely 
the spectral properties of photons emitted in nonlinear 
Compton scattering [28] . See [31] for a discussion of such 
effects beyond the monochromatic approximation. 



fields, tells us to replace external legs with free particle 
wavefunctions. In a background, LSZ reduction trans- 
forms external propagators into incoming and outgoing 
fermion wavefunctions; the following short calculation 
shows this explicitly in plane waves (although the result 
holds more generally). According to LSZ reduction, the 
incoming electron wave function is given by 



¥«(x) = -i /d 4 y S{x,y)[-i$ y -m]e 



(A3) 



Since K $ = 0, we find -i$ -> f + V q ft in (Al). The 
y ± ,y~ integrals set q_ i± = P-,±. Writing q = p + tk we 
get, after simplifying the spin term, 



9?(x) = K px u p e-^- x 
dt 



(A4) 



hi 



2tt 



Vn 



t + ie 



-it((f> x — (f>y)—i J Vp 

4>y 



Now we perform the t-integral. The first term in the 
round brackets gives a delta function which sets <j) y = <p x . 
The second term yields a step function, and the resulting 
(j) y integral is exact. Performing this integral we obtain 
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Appendix A: The Volkov propagator and 
wavef unct ions 



1. LSZ reduction 

We write eC(<p) — a{(j>) from here on to compactify 
notation. We work in the Furry picture, treating the 
coupling to the background exactly and the interactions 
between the quantised fields in perturbation theory as 
normal. Feynman diagrams are therefore built from or- 
dinary QED vertices and the spinor Volkov propagator 
5, which is the inverse ofe-i[0 — ml: 



S(x,y) =i 



d 4 q 



-iq.(x-y) 



k . x 



(2tt) 4 ^f-m + ie 
where we have defined 



(Al) 



2k. p 



V n - ^ , (A2) 



2k. p 



and K = j°K^j°. For ^-matrix elements we also need 
LSZ reduction, which, in the absence of background 



^ n (x) = i^^exp 



<Px 

( — ip.x — i J 



(A5) 



which is the Volkov electron wavefunction (the solution 
to the Dirac equation in a plane wave) with kinetic mo- 
mentum p M and spin a in the far past. Its current is 



1 

2m 



(A6) 



which corresponds to a particle with kinematic momenta 
0. Similarly, see [16], the outgoing electron wavefunc- 
tion which carries kinetic momentum p M and spin a in 
the far future is given by 



^l ut (x) = uZ SK px exp li(p + a 



4>f 

).x-i I SV P ) , (A7) 



in which, and from here on, 5 means a —> a—a^. Positron 
solutions are obtained by sending u —> v and a — >> —a. 
In summary, LSZ transforms external lines into Volkov 
wavefunctions, which describe (on mass-shell) particles 
in a background plane wave. 



2. Normalisation 

The use of the Furry picture Feynman rules and the 
LSZ formulae above correspond to calculating .S'-matrix 
elements using equal-time quantisation as normal, with 
asymptotic momentum states | p ) obeying 



(q\p) = 2p (27r)^ d (p-q; 



(A8) 
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The free states evolve in time to become Volkov wave- 
functions. The structure of the plane wave means firstly 
that these wavefunctions are naturally normalised on 
the lightfront (not as in (A8)), and secondly that S- 
matrix elements of such wavefunctions conserve overall 
p := {p-,p±} : not three- vector p. We give here a clear 
prescription for dealing with normalisations which elim- 
inates the need for volume factors or trying to compare 
the infinite volumes S 3 (p) and 5± _(p), see also [7Q| f7Tj . 

An ^-matrix element calculated using the Volkov solu- 
tions ( A5 ) and ( |A7[ ) for an incoming electron, momentum 
and a set of outgoing particles with momenta {pf} 
takes the form (sums/products of pf are implicit) 



function for the incoming electron: 

S (0) =-i f d 4 x jV+vooUujsiilp - m)^(x) 



= (2nf5^_(p , ^a OQ -p)u : 
f 



X I 



2k + y (A14) 

> 

J d0 (Voo - V)exp (i jv™ - V 



The integral over <\> needs to be regulated. Inserting a 
small convergence factor we find, for a ^ 0, 



5/i = (27r) d ^ )X (p / -p)M(p->p / ), 



(A9) 



S<°> = (2n) 3 5l_(p'+a O0 -p)u p ,^-u p e i 



where 



which defines M. Now, we should really consider scat- 
tering between properly normalised wavepackets rather 
than momentum states. Final states will always be in- 
tegrated out to obtain the full probability with Lorentz 
invariant measure 



0/ 

/< 



{Voo - V) 



(A15) 



(A16) 



E=n 



PS 



f 



d 3 Pf 



(27r) 3 2p /0 



(A10) 



and it is enough to consider only the incoming electron 
wavepacket. This corresponds to multiplying ( A9) by the 
factor 



We now apply (A13). The spin sum and average in our 
case gives 

and it follows that the total probability is 



(A17) 



-/ 



d 3 p' 



(2^2^ 2p_ 



(27r) 3 e_(p'-p)-4p 2 _ 



d 3 p 



with J d 3 P m P )\ 2 = i • (ah) = J ^y-V ff.-tP' - p) = i , 



(A18) 



The S'-matrix element mod-squared then becomes 



1^ 



d 3 p \^j(p)f 



(27r) d (Si |X (p / -p)|M(p-^p / )|- 



(A12) 



2p_ 



Making the usual assumption that the wave packet is 
sharply peaked correspon ds to c alculating ( A9 ) and then 
dropping the first line in (A12) for \Sfi\ 2 . In short, the 
incoming electron should carry a normalisation factor of 
l/2p_, rather than the usual l/2p , at the level of the 
probability, the final expression for which is, summed 
(averaged) over final (initial) polarisation and spin 

F = \ E ^(2 7 r) 3 <5i ) _(p / -p)|M(p^ P/ )| 2 . (A13) 

Pf,a,e 



3. Example: elastic scattering 

As an example, consider the matrix element for elastic 
scattering at tree level. According to the LSZ reduction 
formulae this is obtained by amputating the Volkov wave 



where we used the Lorentz invariance of the measure in 
the second line to change variables. The reason why the 
probability is unity is discussed in Appendix [B] it is a 
manifestation of the IR problem. 



Appendix B: IR structure 

Here we establish in general how soft photons affect a 
given Feynman diagram in the Furry picture. We focus 
on the leading order IR divergences, following Weinberg's 
treatment [33] . 

We use dimensional regularisation to take care of the 
IR divergences, working in 1 + d dimensions with d > 3. 
Noting that our background field singles out a lightlike 
direction (j) = k.x ~ x + , we follow [72J and place the ex- 
tra dimensions into the d—1 directions transverse to the 
background. This means in particular that all the struc- 
ture of the plane wave is preserved by the regularisation. 
The measures in position and momentum space are then, 
in lightfront coordinates, 



dx :- 
dq:~- 



dx+dard'*- 1 ^ = 



d d+i 



2k + 
dq + dg_dgf _1 



dx-d^x 



(Bl) 



(27T) 



d+1 



7T 



(27T) 



9 



As we are only interested in the IR, we assume there is a 
cutoff in place to take care of the UV divergences, regard- 
ing which we note the following. By employing a proper 
time representation, all propagators can be expressed in 
terms of heat-kernels, and these are easily continued to 



d ^ 3 following [73]. The short-time expansion of the 
heat-kernel then gives a very convenient method for iden- 
tifying UV divergences even when background fields are 
present. See [63] for an example. 



1. Soft photon correction to an external line 

We concentrate on incoming electron lines, other external lines can be treated similarly. A Feynman diagram in 
which the incoming electron emits n soft photons (which may be real or virtual, we consider both cases below) with 
small momenta lj contains the term 

E n (x n+1 ) :=(-ie) n J dx n ..A Xl G(x n+1 ,x n )^e a ^ . . . G(x 2 , x^e^V^x^ . (B2) 

We will first factor out the infrared divergent part. To begin, consider n = 1, 



2 1 

E 1 (x 2 ) = e Jdx 1 Jdq K q2 g2 ^ ^ 2 + .^ K ql ^K pl u p exp ^ - iq.x 2 + i(q + h -p)-Xi - i J V q - i J V p 



(B3) 



The x~ and x^-integrals sets q = p — 1 1 . (In the absence of the field, we would also be able to perform the x + integral 
to set = p^ — .) Since we are only interested in the soft sector, and in particular the divergent terms, we employ 
the usual eikonal approximation, replacing k.q — >> k.p in K. Define now a new variable t by q 2 — m 2 = 2k. qt « 2k.pt. 
Changing integration variable from q + to t we find, again to lowest order in I 

2 2 

= 2^p / ^ / ~^ Kp2 ^ + ie + ^ KpiJ^KpiUpexp ^ - it(<j> 2 - <t>i) - i(p - h)-x 2 - i J V p - i J . 

1 

(B4) 

We have condensed our notation further: when tt appears under a (j) integral it means 7r(p;0), and when it appears 
outside such integrals it means 7r(p; 00). Consider the term in square brackets. Performing the t integral, the ^-term 
leads to a delta function setting (j) 2 = 0i. This is the contribution from the lightfront zero mode [61 , which is 
interesting in itself but does not contribute to the IR divergence and so we drop it. The remaining term in the square 
brackets gives 9((j) 2 — 0i), which restricts the <j)\ integral. The spin term can be simplified to (^-\-m)Kj^Ku = 2n^u 1 
and we then find 

02 02 

E 1 (x 2 ) = *™(x 2 ) e il ^ j exp ( - i J ^) . (B5) 

— OO 0! 

The calculation is easily extended, so that we obtain 

n x 2 n x 

E n (x) = % n (x) exp (i lj-*) ("^) / d^n • • • J dfa f[ Trf exp ( - i J . (B6) 



The leading term is the incoming Volkov solution. The evaluation of the remaining integrals, which contains the 
IR divergence, depends on whether a) the leg is attached to a 'hard' vertex and is therefore part of a multi-particle 
scattering process or b) the leg continues to an outgoing line, in which case there are only soft photons in the process. 
We now consider these two cases, which are illustrated in Fig. [4] 



2. Between a soft and a hard place 

We begin with the left hand diagram in Fig. [4j considering the effect of the soft photon lines in (B6) (which can 
correspond real emission or virtual loops, as we will need both) on a hard scattering process. We will see that in this 
case, the IR divergence does not depend on the properties of our background field. 
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FIG. 4. Left: an external leg emits soft photons (or emits and absorbs virtual photons) as part of a scattering process with 
a hard vertex at Right: emission of many soft photons from a single electron. There is no hard scattering part to this 
diagram. 

So, we assume that E n (x) is connected to a hard vertex. In that case, there is a hard photon momentum in the 
exponent of (B6) at position x, and relative to this we can neglect the term lj. Since we have reduced the spin terms 



to products of 7r's in (B6), the integrand therein is symmetric with respect to the except for the step functions. 
This simplifies when we sum over all the orders in which the soft photons can be emitted (i.e. when we account 
properly for all diagrams): 

n x I x \ 

£ E n { X ) = (~YU [ exp -i / ^ ), (B7) 

l-perm V * ' J = l_ oc \ j I 

and we see that the step functions drop out. The divergent part of the ^-integrals comes from the region before 
the pulse turns on. Since the leg is attached to the hard vertex, the upper limit of the integrals in the exponent is 
unimportant as long as it is finite, and we can make the following replacement without affecting the leading divergence, 

x 

f h - ^ j <f>j ■ (B8) 

-co -co 

Performing the <f> integrals with the help of a convergence factor as usual, we find 

E^=wn^-- ( B9 ) 

j=l 

in which the soft contributions factor off, meaning that the Feynman amplitude factorises into a hard part and a soft 
correction, as for QED without background fields. Similar expressions hold for other external lines. Hence, we have 
reduced the IR problem in plane waves to the case with no background, with no surprises thus far, and we may proceed 
as in [33, §13] to show that the leading infrared divergences cancel to all orders when one sums the probabilities for 
indistinguishable processes. These statements hold for processes in which a 'hard part' can be identified, i.e. assuming 
that the external lines are connected to hard vertices. We turn now to single electron processes with only soft vertices. 

3. Soft processes 

The ^-matrix element for a diagram with a single electron and n soft vertices, see the right hand diagram of Fig. [4j 
contains 

S n := -ie J dx ^{x)^e iLx E^x) . (BIO) 

Since the photons are soft, k.q ~ k.p and k.p' ~ k.p. The spin terms can then be simplified and taken outside the 
integrals, leaving 

co n 2 n j 

S n = (2n) d 5 3 _ t± (p' + aco - p y e u p ,^u p (-^j J d0„ J 64> n -! . . . J tyi f[ < j exp (i J , (Bll) 

-co —co —co 3~ 1 
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with as in (A16). (We drop the factor Y^h fr° m the delta functions. In a more rigorous treatment the soft photon 
energies should be restricted so that this sum is less then some given energy, see [33 . The divergent part is still the 
same.) As before, S n simplifies when we sum over the permutations of the lj. 



]T S n = (2n) d 5 3 _J P > 



I— perm 



9 _/ It 

u ——u 
2k + 



le 
kp 



n 



TTj exp 



"3 "3 




(B12) 



We have derived this formula for n > 1 (with Eq = \I/ m ), but it also holds for n = 0, when it describes elastic scattering 
at tree level, see Appendix | A 3| above. The significant difference compared to the case of hard-soft factorisation is that 
in soft processes, the outgoing electron's momentum is fixed by classical momentum conservation, in other words by 
the properties of the background field, and in particular a^. 



The soft photons can be real or virtual. For each real emission we multiply (B12) by a polarisation vector e, giving 



le 
k.p 



e.7rexp 




ee 



7T 



P_ 
Lp 



(B13) 



At the level of the probability we sum over polarisa- 
tions, which gives minus the above expression squared, 
and then integrate over the photon momenta. We get 
the same factor for each photon, but with a symmetry 
factor of l/n\. The contribution from the emission of n 
soft photons is therefore 



nl 



7 



d d l 1 

(2n) d 2f \l.7T 



P_ 
Lp 



For each virtu al photon we choose lj 
multiply QB12D by 



(B14) 



-L = / and 



/ 



dl 



I 2 



(B15) 



We then have 



/ 



dl- 



P + ie 



k.p 



^TTiTTj exp 



> 



and 



< 



(B16) 
and 



We divide this into two parts 
change variable I —I in the <-part. We can t hen close 
the /o-contour in the lower half plane, and (B16) becomes 



d d l 1 



e 

k.p 



?<e + ) , 
(B17) 



with obvious notation. The imaginary part of (B16) di- 
verges like 



Im 



/ 



d<t>, 



(B18) 



but we will shortly see that it drops out of probabilities. 
In the real part, the fa integrals can be performed and 
are finite, and (B16) becomes 



d d l 1 



p_ 
Lp 



(B19) 



(27r) d 2/ 

We get one such factor for each virtual photon with a 
factor of l/2 n n!, which is the number of identical permu- 
tations of the sum over lj for n virtual photons. Summing 
over all n we get the all-orders loop contribution to a soft 
process, 



exp 



d d l 1 

(2ir) d 2f \l.7T 



P_ 
Lp 



-hi... 



(B20) 



This contribution appears mod-squared at the level of 
the probability, which removes both the leading factor 
of one half, and the divergent imaginary term; the lat- 
ter are the usual phase divergences. Hence, returning to 
the example of Appendix |A 3[ the probability for elastic 
scattering including all soft loop contributions is given 
by the modulus squared of ( B20| ), 



exp 



d d l 1 

(2ir) d 2f \l.7T Lp 



(B21) 



When ^ 0, i.e. when the background field is unipolar, 
the loops give an IR divergent contribution as d — >• 3. 
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